Seeing Anderson Localization 



Dpt. of Physics, 



M. Hilke 

McGill University, Montreal, 



Canada H3A 2T8 



Anderson localization was discovered 50 years ago to describe the propagation of electrons in 
the presence of disorder The main prediction back then, was the existence of disorder induced 
locaUzed states, which do not conduct electricity. Many years later it turns out, that the concept of 
Anderson localization is much more general and applies to almost any type of propagation in time or 
space, when more than one parameter is relevant (like phase and amplitude). Here we propose a new 
optical scheme to literally see Anderson localization by varying the optical wavelength or angle of 
incidence to tune between localized and delocalized states. The occurrence of Anderson localization 
in the propagation of light, in particular, has become the focus of tremendous interest due to the 
emergence of new optical technologies and media, such as low dimensional and disordered optical 
lattices Q y] . While several experiments have reported the measurement of Anderson localization of 
light 0,lli,i, many of the observations remain controversial because the effects of absorption 
and localization have a similar signature, i.e., exponential decrease of the transmission with the 
system size 0, [!]• In this work, we discuss a system, where we can clearly differentiate between 
absorption and locaUzation effects because this system is equivalent to a perfect filter, only in the 
absence of any absorption. Indeed, only one wavelength is perfectly transmitted and all others 
are fully localized. These results were obtained by developing a new theoretical framework for the 
average optical transmission through disordered media. 



In order to literally see Anderson localization, we con- 
sider the system shown in figure 1. The setup is com- 
posed of Nf filters in series. Each filter is composed of 
Ni random optical layers and each layer has a refractive 
index rij and a thickness dj as defined in figure 2. For 
optimal results, it is important that these layers are very 
well defined. This is possible, for instance, by using high 
accuracy multilayer growth techniques, such as molecu- 
lar beam epitaxy (MBE). Good material choices include 
large band gap materials like GaN (A = 3.2eV, n :^ 2.3) 
and InN (A = 2eV, n 3.1), which can be grown by 
MBE with atomic precision [13] • The absence of any 
surface roughness ensures that this multilayered system 
is one-dimensional in nature. Surface roughness would 
induce scattering at stray angles, which leads to a ID 
to 3D dimensional crossover , thereby strongly reduc- 
ing localization effects, since localization in 3D is much 
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12| and therefore much more chal- 



weaker than in ID _ 
lenging to observe [5]" 

To illustrate the general behavior of our system, we 
considered the simplest random system, a binary dis- 
tribution, where we assume two kinds of materials, one 
with a refraction index of no = 1 and the other with 
ni = 3. Adding more materials will not change the re- 
sults qualitatively and the analytical expressions derived 
below remain valid for any choice or combination of mate- 
rials. However, to maximize the effect of disorder (small- 
est localization length) the variance of refraction indices 
should be largest. Since in most materials the refraction 
index rarely exceeds n = 3, an efficient way to increase 
the variance is to consider a discrete distribution. Indeed, 
for a binary random distribution {no = 1 and ni = 3} 
the variance is one as compared to a uniform random 
distribution 1 < n < 3, where the variance is reduced to 
one-third. 
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FIG. 1: The top shows the arrangement of the filters, where 
a grating can be inserted in order to measure the wavelength 
dependence of the transmission. This transmission is shown 
in the bottom left figure when the grating is placed after 1, 3, 
10, 30, or 100 filters, which increases the quality factor filter 
with the number of layers. The bottom right figure represents 
the dependence of the transmission on the number of layers 
at several wavelengths (480, 497, 499, and 500)nm, indicating 
an exponential decay 



The filters are placed in a way to allow the insertion of 
a grating between two of them. The grating is necessary 
to resolve the spectral composition of the transmitted 
light. The position of the grating can be changed in order 
to measure the transmission after 1, 2,. . . , or Nf filters 
in series. When using a white light source, this allows 
to spectrally resolve the transmission after any number 
of filters. Anderson localization will then lead to an ex- 
ponential decay of the transmission with the number of 
filters (see figure 1). 
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FIG. 2: The top picture illustrates the details of a random 
sequence of layers composing one filter. In the bottom left 
figure the transmission as a function of wavelength is shown 
for a system of 200 layers averaged 10^ times. The analyti- 
cal result from equation ^ is shown with a continuous line. 
The numerical average and relative fluctuations are shown 
with square and triangles. The transmission as a function 
of wavelength for different incidence angles (0, 15 and 30°) 
is depicted in the bottom right figure, where the continuous 
lines correspond to the analytical expression JS}, whereas the 
dotted lines correspond to the numerical log-average. 



A major difficulty in observing Anderson localization is 
the existence of strong fluctuations In general, the 

relative fluctuations of the transmission {ST/ (T)) diverge 
with decreasing transmission, where 5T = ^J {T — {T)Y 
and (T) is the average transmission. This is a fundamen- 
tal problem, which can be circumvented by considering 
ln(T) (the transmission in decibel) or the inverse trans- 
mission since in both cases the relative fluctuations de- 
crease with decreasing transmission as shown in figure 2 
for ln(T). Therefore, in order to reduce the fluctuations 
we can average the logarithm of the transmission over 
several disorder configurations, which converges much 
faster than an average over the transmission. The nu- 
merical examples shown in figure 1 are based on a small 
log- average over 100 configurations and assuming 100 dif- 
ferent filters, each with an average of 50 random layers 
and a spectral resolution of Inm. Such a log-average 
over a small number of configurations is feasible exper- 
imentally, since it simply involves averaging the signal 
in decibel over a few configurations. In comparison, a 
higher averaging is used in figure 2, where 10^ configura- 
tions with 200 random layers was used, leading to even 
less fluctuations. Considering a system of Nf different 
filters allows one to perform a configurational average 
over the disorder, where the space between two filters 
is equivalent to an additional layer, as long as the dis- 
tance between two filters follows the resonance condition 
described below. However, this is equivalent to simply 
considering a single filter with more random layers and 
then performing a configurational average, which is what 



we use in the remainder of this letter. 

The results were obtained using a standard method 
to describe a multilayered optical system with a transfer 
matrix. Each matrix describes the transmission after one 
layer, where we assume the layers to be normal to the 
X direction. For TM (transverse magnetic) waves, where 
the magnetic field component is parallel to each layer 
and the electric field component is at an angle (the 
angle of incidence) to the first layer, the transfer matrix 
can be written as [3]: 



cos{(j)j ) —i sin{(j)j ) / 
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Mi 



(1) 

For an incoming wave, of wavelength A, the field com- 
ponents on one side of a stack of Ni layers can therefore 
be related to the fields on the other side by taking the 
product of the transfer matrices (M = Yl^liMj). As- 
suming a vacuum impedance zq before and after the filter, 
the transmission of light is given by T = 4|Mii +7M12 + 
^21/7+^22 1 where M^i are the matrix elements of M 
. The matrix elements of Mj used in equation ([1]) are 

material dependent with = (27idj / X)^ n^^ — sin^(6>) 
and7_^- = (zo/n^)yri? 



sin (6) . For TE (transverse elec- 
trical) waves these expressions have to be replaced by 
Ey -Ry, and 7^- = ZQ^^n]-sm^{e). 
Absorption exists when the imaginary part of the re- 
fraction index is non-zero (^{n^} 7^ 0). ^{^j} can 
be substantial when photon energies approach the band 
gap energy or when the materials are conducting. With 
absorption the localization-delocalization transition, de- 
scribed below, cannot be observed. For instance, for 100 
layers and an imaginary part of the refraction index of 
only 0.01 the transmission at the resonance is decreased 
by about a factor 100 so that absorption effects can un- 
ambiguously be distinguished from localization effects. 
This is crucial, and in the remainder of this letter we 
will only consider the case without absorption, for which 
the localization-delocalization transition can be obtained 
under the following condition: 

For a given layer there exists a set of thicknesses 
dj = mXo/2nj (m any integer) for which the layer is 
totally transparent at normal incidence and for the res- 
onant wavelength Aq. This corresponds to sin(^j) = 
and to Mj equal to the identity matrix. When randomly 
combining different layers with the same resonance con- 
dition, the overall resonance condition is preserved. The 
position of these resonances was discussed by Cristanti 
[15]. Away from the resonant wavelength the system be- 
haves like a random system, leading to Anderson local- 
ization. This type of disorder was also studied in the 
context of electronic transport in ID and 2D 1^, 17, lit . 
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where metal-insulator type transitions were found. Fig- 
ures 1 and 2 show this resonance, where the transmission 
is one (transparent) at Aq = 500nm and then decreases 
exponentially (T e~^''^~^^^ ) away from the resonance 
condition. K is a form factor which is proportional to the 
number of random layers and also depends on the vari- 
ance of the disorder as discussed below. This dependence 
on the number of layers can then be used to tune the qual- 
ity factor of the filter. Without this special resonance 
condition, the transmission through random optical lay- 
ers decays exponentially at all wavelengths [6l, |l9|, l20| . 

We now turn to characterize the transmission through 
this random system by evaluating the product of M^'s 
composed of random elements. Similar products of ran- 



dom matrices have been widely studied and lead to ma- 
trix elements, which increase exponentially with the num- 
ber of products. The rate of this exponential dependence 
is termed the Lyapunov exponent and will be derived 
analytically below. The Lyapunov exponent also corre- 
sponds to the inverse of the localization length, beyond 
which the transmission vanishes. To obtain analytical 
results the trick is to consider the square of the fields, 
\Hj\'^ instead of Hj. This allows to concentrate on the 
dominant behavior beyond the plane wave oscillations. 
A similar technique was actually used in the context of 
electronic transport [2l|, [22|, 113] and we extend it here to 
optical systems. Squaring equation ([T]) then leads to the 
following iterative equation: 



\Ey_,\' = (sin(</.,)7,)2 cosHcP,) -icos{cP,)sm{^,h, ] { I^JP . (2) 

V ^ V V ^ 



For Ni layers, the total transmission is now deter- 
mined by the product Fq = (U^li Mf''^)FNi, which 
can be averaged over the disorder, leading to, (Fq) = 
(M^><^)^^(F^), where (M^><^) is the disorder average of 
Mj^^ assuming that the material parameters of the lay- 
ers are not correlated. The disorder average depends on 
the distribution of the parameters entering Mj^^. The 
leading behavior is obtained by taking the eigenvalues of 
^M^^^), from which the Lyapunov exponent can be ex- 
tracted as A = max{ln |Eig((M^><^)|)}. This yields an 
expression for the average resistance (inverse transmis- 



sion), where (1/T) 



Hence, perfect transmission 



corresponds to the case where A = 0. Interestingly, this 
expression can also be related to the log-average trans- 
mission (InT) = —ANi/2 and to the average transmission 
(T) = e~^^'/^, where the factors of 1/2 stem from the 
properties of the distribution of T's [22|. The next step 
is obtaining (M^^^) for a given distribution. For discrete 
distributions, like the binary one this is quite straightfor- 
ward, since the layer properties are determined by only 
two possible configurations, {do, ^0} and {di, ni}, re- 
spectively. This yields 



(M^><^) = (Mo^><^ + Mf><^)/2 



(3) 



where Mq and M^^^ are simply M^^^ 



with {dj = 

do^ rij = no } and {dj = <ii, = ni }, respectively. 
For continuous distributions, analytical expressions for 
(M^^^) are often too complicated to be useful and can 
be integrated numerically instead. 



Comparing these analytical expressions to the numer- 
ical results demonstrate a remarkable agreement as seen 
in figures 1 and 2. The analytical expression using equa- 
tion (|3]) reproduces all the main features and shows the 
resonance behavior at A = Ao- Evaluating the logarithm 
of the eigenvalues of expression ([3]) close to the resonance 
gives A ~ (A — Ao)^, which is perfectly reproduced by 
our numerical results and determines the form factor 
which depends on the variance of the disorder. (Zero vari- 
ance leads to = or A = 0). Deviations exist when 
looking at the angle dependence, as shown in figure 2. 
While the analytical expression predicts two transmission 
peaks with a weak minimum in between, the numerical 
results show that the minimum is more pronounced. This 
can be attributed to restricting the average to the second 
moment, which can lead to differences with the numeri- 
cal results and in some cases even lead to fiuctuations in 
the Lyapunov exponent [23] . 

Quite remarkably, it is now possible to tune the wave- 
length of the resonance, simply by tilting the filter. We 
present the angle and wavelength dependence of the 
transmission in figure 3 for both the TE and TM waves 
using 200 random layers averaged over 100 configura- 
tions. Two branches can be observed, one which is per- 
fectly transmitting and a second one which vanishes with 
increased angle. In addition, there is an angle 6 :^ ±1.25 
with perfect transmission for the TM incidence, where 
O = arctan(ni/no) rad and ni/no = 3 for our particular 
configuration. This angle corresponds to the Brewster 
angle. 
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FIG. 3: The wavelength and incidence angle dependence of 
the transmission through a random multilayer on a linear 
color scale for TE waves (left) and TM waves (right). 
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FIG. 4: The transmission is shown on a linear color scale as 
a function of the angle of incidence in the y and z directions. 
Each subgraph spans -7v/2 to +7r/2 and corresponds to a dif- 
ferent wavelength. The resonant wavelength (Ao = 500nm) 
corresponds to the bottom left panel. The horizontal (verti- 
cal) axis would correspond to the angle dependence of a TM 
(TE) wave. 



By combining the expression for TE and TM waves it 
is possible to obtain the results for any polarization, sim- 
ply by using a linear combination of TE and TM waves. 
Moreover, we can analyze what happens when tilting the 
filter in one or the other direction. The tilt angle can 
be varied in two spatial directions, one corresponding to 
the TM direction and the other to the TE direction, or 
a combination of both, which is seen in figure 4. Per- 
fect transmission is now identified in a form of a ring, 
whose diameter depends on the wavelength. This angu- 
lar dependence therefore allows for a direct and beautiful 
visualization of the localization-delocalization transition 



and hence of localization. 

Summarizing, we have shown a new filter design, where 
we can directly visualize Anderson localization by using 
a localization-delocalization transition. The filter design 
is perfect in the sense that only one visible wavelength 
is transmitted perfectly and all others are exponentially 
suppressed. Not only does this provide for an optimal fil- 
ter design but also allows us to see Anderson localization 
50 years after its discovery. 
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